Minimal energy packings of nearly flexible polymers 
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We extend recent studies of the minimal energy packings of short flexible polymers with hard- 
core-like repulsions and short-range attractions to include bond-angle interactions with the aim of 
describing the collapsed conformations of 'colloidal' polymers. We find that flexible tangent sticky- 
hard-sphere (t-SHS) packings provide a useful perturbative basis for analyzing polymer packings 
with nonzero bending stiffness only for small ratios of the stiffnesses for the bond-angle and pair in- 
teractions, kl^^^ /kc < 0.01 for A/" < 10, and the critical ratio decreases with the number of monomers 
N . Thus we argue that future studies should not use flexible t-SHS minimal energy packings as 
a reference set to enumerate and characterize the structural and mechanical properties of polymer 
packings with finite bending stiffness. 



PACS numbers: 64.70.km,81.16.Dn,02.10.Ox,87.15.hp 



Protein folding and other examples of polymer col- 
lapse in dilute solutions are complex processes that in- 
volve the cooperative motion of thousands of atoms. A 
number of early studies of polymer and protein struc- 
ture and dynamics employed hard-sphere models with 
only steric interactions, bond-length, and bond-angle 
constraints to understand polymer elasticity [1] as well 
as secondary structure [2] and packing in hydrophobic 
cavities in proteins. More recent studies have imple- 
mented Monte Carlo, Brownian dynamics, and molecular 
dynamics (MD) simulations that include excluded vol- 
ume and stereochemical constraints as well as solvent- 
mediated attractive interactions to investigate numerous 
examples of polymer collapse. However, recent advances 
in colloidal synthesis and self assembly have enabled stud- 
ies of "colloidal polymers" with hard-core-like repulsions, 
short-range attractions, and greater chain flexibility than 
typical synthetic or biological polymers ^S^, which has 
prompted renewed interest in simple polymer models. 
We seek to develop the ability to theoretically predict 
which structures are most likely to form for a given set 
of interactions and solvent conditions, which will reduce 
the need for trial-and-error synthesis of self-assembled 
nanostructures. 

The ground state packings of short, flexible tangent 
sticky hard sphere (t-SHS) polymers have recently been 
characterized via complete enumeration techniques. jH 
These packings possess a wider range of symmetries and 
shapes than the ground state packings of polymers with 
longer-ranged attractive {e.g. Lennard-Jones|7|) pair po- 
tentials. Further, the number A^micro of distinguish- 
able, energetically degenerate ground-state flexible t- 
SHS packings ( "microstates" ) grows exponentially with 
the polymerization index A/", with a corresponding in- 
crease in the diversity of paths followed by the cova- 



lent backbones. In this manuscript, we will character- 
ize the ground state packings for short, semi- flexible t- 
SHS polymers with finite bending stiffness /c^ using the 
ground-state flexible t-SHS polymer packings as a refer- 
ence set. In particular, we will calculate the potential 
energy of semi-flexible polymer packings as a function 
of kb and identify the maximum bending stiffness k^^^^ 
below which the reference set of t-SHS packings are the 
lowest energy structures. These studies will inform fu- 
ture exact-enumeration studies of polymer packings with 
larger bending stiffness. 
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FIG. 1: Schematic of the low- lying minima in the energy land- 
scape for polymer packings as a function of the configurational 
degrees of freedom {r}. For vanishing bending stiffness kb = 
(left), the ground-state flexible polymer packings (stars) are 
energetically degenerate and correspond to the lowest Amicro 
inherent structures (minima of the black curve). Increasing kb 
(right) deforms the energy landscape and raises the energies 
of the Amicro minima that correspond to the flexible poly- 
mer packings. For kb < k^^^^ , the Amicro lowest energy semi- 
flexible polymer packings correspond to the Amicro ground- 
state flexible polymer packings. 

Figure [1] summarizes our theoretical approach. Con- 
sider a polymer model with harmonic pair and bond- 
angle potentials with respective "spring constants" kc 
and and equilibrium bond angle Oeq. The ground- 
state flexible {kb = 0) t-SHS polymer packings can pro- 
vide a useful reference set of configurations for studying 
the lowest energy semi-flexible {k^ > 0) t-SHS polymer 
packings when there is a one-to-one mapping from the 
A^micro ground state polymer packings for kb = to the 
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^micro lowest energy polymer packings with finite bend- 
ing stiffness. More specifically, our approach assumes 
that the A^micro zeio-kb t-SHS packings lie within the 
basins [8] of the A/'micro lowest-energy packings at finite 
as shown in Fig. [T] (right). Increasing kb > deforms the 
energy landscape, raises the energies of the A^micro flexi- 
ble polymer packings, and breaks their degeneracy. For 
kb < ^5^^*, the A/'micro lowest energy semi- flexible polymer 
packings correspond to the A^micro ground-state flexible 
polymer packings. In contrast, for k^, > k^^^^ the energy 
landscape is sufficiently deformed that new low-energy 
packings (that do not belong to the set of k^ = ground- 
state polymer packings) are among the N micro lowest. 
For example, in Fig. [T] with A'micro = 8, if the energy 
of polymer packing 9 drops below any of the energies of 
packings 1 through 8, kh > k^^^^. Below, we will calcu- 
late the total energy (that includes bending and attrac- 
tive contact energy terms) for the ground state t-SHS 
polymer packings as a function of k^ and determine k^^^^. 

The pair potential for t-SHS^ polymers with spherical 
monomers of diameter islj 



oo ^ r < D 

-e , r = D 

^ r > D ^ noncovalent 

oo , r > D , covalent. 



(1) 



where r is the center-to-center separation between 
monomers and — e is the contact energy. In Eq. (TJ 
the term "covalent" ("noncovalent") refers to monomers 
that are (not) chemically adjacent and permanently con- 
nected. All A'micro distinguishable polymeric paths (mi- 
crostates) through A^-monomer, A'c-contact t-SHS pack- 
ings are obtained via the complete enumeration proce- 
dure described in Ref. [6]. Values of A/'micro for ground- 
state (maximal- A^c) t-SHS packings of 6 < A^ < 9 
monomers are given in Table HI 



TABLE I: Number of microstates A^micro for maximally con- 
tacting t-SHS packings, where A^micro for A" > 6 includes 
a factor of two associated with packings possessing chiral 
enantiomers. [10] 



A^ 

6 

7 



Nc 
12 
15 
18 
21 



cro 

50 
700 
6429 
122060 



For numerical calculations, we employ a continuous 
version of /7ss(r) with short-range attractive interac- 
tions and harmonic bond-length constraints illustrated 
in Fig. El (a). 



kc / r \ ^ 
e+y(^-l) , r<r. 



Ue{r) 



(2) 



where Tc/D = oo for covalently bonded monomers and 
Tc = = D{1 -\- ^y2e/kc) for noncovalently bonded 
monomers. Note that /7ss(r) is the kc ^ oo limit of Uc 
in Eq. O To model finite bending stiffness, we employ 
a harmonic bond-angle potential used in many previ- 
ous computational studies of organic molecules, peptides, 
and proteins [Tl|, [l2| : 



TP 

^bend 



(3) 



The bond angle between adjacent monomers (z, z + 1, 
i-\-2) is defined as = tt — acos[(6i • 6^+i)/| 6^6^+1 1], with 
bi = fi+i — fi. /7^end minimized at the equilibrium 
bond angle 0'^ = 0^^^ and k^ is expressed in units of e. 
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FIG. 2: Schematic of (a) the short-range attractive and har- 
monic bond-length potentials in Eq. [2] and (b) definition of 
the bond angle for the harmonic bond-angle potential in Eq. 

ini 

The total potential energy of an A^-monomer chain is 



ch am 



bend* 



(4) 



For /c5 = and the limit /Cc ^ oo, configurations corre- 
sponding to local minima of t/chain are clearly identical 
to those obtained in studies of t-SHS packings with no 
bond-angle interactions, i.e. all A/"- monomer, A/'c-contact 
microstates have /7chain = —NcC. Bond-angle interac- 
tions break this degeneracy; t/chain increases by different 
amounts for each microstate. 

In the sticky hard sphere {kc oo) limit, to first or- 
der, increasing the bending stiffness does not change the 
structure of the polymer packings; in the limit k^/kc 
0, the increase in energy Uchain{kb) - Uchain{kb = 0) 
scales linearly with kb. Each reference {kb = 0) mi- 
crostate possesses a set of A/" — 2 "reference" bond angles 
{l9[''^, 6>2''^, e';^^^}. Thus, the leading order (LO) 
prediction for the energy of the jth microstate is 



^chain,LO 



where 



N-2 



kb 

2^2^ 



(5) 



(6) 



and we assume $1^ = O^q for all i. Since all cj are 
available from the flexible t-SHS packings, LO predic- 
tions for the energy of polymer packings for arbitrary kb 
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and O^q can be made using only information from the 
flexible reference polymer packings. Predictions for the 
minimum, maximum, and average energies of polymer 
packings {Umin, Umax, Uavg) from Eq. [5] for N = S and 
two physically relevant Oeq are shown in Fig. [3l Alkanes 
and other polymers with tetrahedral bond-angle order 
have Oeq ^ 9tet = acos(-l/3) = 109.47°, while Oeq = tt 
corresponds to alkynes and other polymers with linear 
bond-angle order. The bond-angle energies are higher 
for Oeq = TT siucc the compact structure of small- t- 
SHS packings allows few "straight trimers" (angles with 
0' = 7r). 




packings as shown in Fig. IH The t-SHS packings possess 
a wide range of bond angles, leadin g to large A^, and 
differently structured packings [si [iqI,1i3| and polymeric 
paths through these packings lead to large variations in 
(0). The variation in both (0) and AO increases rapidly 
with N. 




kb/kc 10- 



FIG. 4: Bond-angle dispersion for t-SHS polymer packings. 
Each data point represents the mean bond angle and the 
rms deviation for a single t-SHS polymer packing. Large 
black, medium blue, and small green circles represent data for 
A = 6. 7, and 8, respectively. 



FIG. 3: Comparison of the leading-order expressions for the 
energy [/chain of A = 8 semiflexible polymer packings in Eq.[5] 
(dashed, Umin] solid, t/avg; dotted, t/max) to the numerical re- 
sults for the energy of minimized t-SHS packings in the pres- 
ence of bond- angle interactions (downward triangles, Umin'-, 
circles, Uavg; upward triangles, t/max) as a function of kb/kc. 
Results for Oeq = Otet and tt are shown in green and blue, 
respectively, t/chain/e = —{Nc — 1) and [/chain = are respec- 
tively indicated by the black and red horizontal lines. 

To validate the leading-order expressions for the en- 
ergy of semiflexible polymer packings in Eq. [5l we em- 
ployed energy minimization techniques (described in the 
Appendix) to relax the t-SHS configurations in the pres- 
ence of bond-angle interactions. The leading-order ex- 
pressions for the minimum f/min, average /7avg, and max- 
imum [/max chain energies closely agree with the nu- 
merical results over a wide range of k^/kc as shown in 
Fig. [3l Both the numerical results and leading-order ex- 
pressions for [/chain are nearly linear in kij/kc over the 
full range studied. We know that k^/kc > k^^^^/kc for 
^chain > siucc frecly-rotatiug chain conformations with 
no pair contacts and 0'^ = Oeq for all bond angles have 
^chain = 0. Thus it is clcar that the critical bending 
stiffness k^^^^ (above which the A^micro lowest energy semi- 
flexible polymer packings are not the same as the A'micro 
ground-state t-SHS packings) is below k^^^^/kc ~ lO""*^'^ 
where [/chain > -NcC. 

For A^ = 8, we find that k^ exceeds /cg^*^ at a rela- 
tively small value for the bending stiffness. Note that 
for l9eg = Otet, (t^max + Ncc) / {Umin + A^^e) 7 bccausc 
there is a large dispersion of bond angles in the t-SHS 



A more refined estimation of k^^^^ can be obtained as 
follows. The A'micro (A') flexible t-SHS packings all pos- 
sess Nc{N) pair contacts (Table |T]) and energy —Nc{N)e. 
Any polymer packing with fewer than Nc{N) pair con- 
tacts must have [/chain ^ — (A'c — 1)^ since the bending 
energy is strictly positive. Conversely, any t-SHS poly- 
mer packing with A'c pair contacts must correspond (in 
the sense of Fig. [1]) to one of the A^micro(A^) flexible t-SHS 
packings. Therefore kl^^^{N^Oeq) is set by the largest 
Cj{N^Oeq) (Eq. [6]) for these packings and the condition 
that the reduced bending energy is less than unity, and 
thus 

kf\N, 0,,) < 27rVc7'^"(iV, Oe,). (7) 

As shown in Fig. [5l k^^^^ decreases monotonically with 
increasing A^ for all Oeq > 7r/3, which stems from the 
increasing angular dispersion with A" shown in Fig. HI 
In addition, for increasing Oeq^ k^^^^ reaches a peak at 
^max that increases with A", and then decreases steadily 
for Oeq > 6>max- For the A^ considered here, the peak 
occurs at a rather small ^max ^ 7r/2 because of the com- 
pactness of the flexible t-SHS reference polymer packings. 
Nonetheless, we expect the key result of decreasing /cg^*^ 
with increasing A^ to continue for A^ > 9. 

To check the analytic prediction of Eq. [71 for /c^^^*, we 
performed molecular dynamics simulations of semiflex- 
ible t-SHS chain collapse for A^ = 8. Monomers in- 
teract via the pair and bond- angle potentials in Eqs. 
[2] and [3l with kc = 3200e and a range of Oeq and 
k}). For each Oeq^ k^, and kc, systems of A's = 1000 
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ble polymer packings with kij > /cg^** using Monte Carlo 
or molecular dynamics simulations with advanced sam- 
pling techniques. 
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'{N, Oeq) from Eq.Elfor AT = 6 (black 
line), 7 (blue line), 8 (green line), and 9 (red line) and MD 
simulations of polymer collapse for N — S (green circles). 



independent self-avoiding random walk polymer config- 
urations were slowly quenched from high temperature 
{ksTi = 4e) to ksTf < 10~^^e using an isokinetic ther- 
mostat [14] and quench protocol T = exp [— t/(10^r)], 
where r = ^/mD'^ /e is the unit of time and m is the 
monomer mass. An overestimate of k^^^^ can be obtained 
by measuring the lowest /c^ at which the minimum energy 
collapsed states obtained via MD have chain energies be- 
low t/min, which cannot correspond to any of the A^micro 
flexible t-SHS inherent structures. 

A comparison of k^^^^ from Eq. [7] and estimates from 
the MD simulations are shown in Fig. [5] for several Oeq. 
In all cases, the appearance of polymer packings with 
energies lower than Umin coincides with contact-breaking 
{i.e. the minimum energy semiflexible polymer packings 
from MD simulations possess fewer pair contacts than 
the flexible t-SHS inherent structures) and occurs for k^ 
that satisfy [UmaAK) + (3^ - 6)] (2 - 3)e. Thus, 
energetically favorable rearrangements of t-SHS polymer 
packings due to finite bending stiffness require breaking 
two or three pair contacts. 

We conclude that under our very strict condition that 
all of the A^micro lowest-energy zero-A:^ t-SHS states 
must maintain a one-to-one correspondence to the A/'micro 
lowest-energy states for finite /C5, the flexible t-SHS min- 
imal energy packings are of limited utility as a basis for 
perturbative studies of polymers with finite bending stiff- 
ness. A less strict condition, e.g. requiring only a one-to- 
one correspondence for ground states, would allow exten- 
sion of perturbative approaches to larger k^. Nonetheless, 
our results suggest that even very small bending stiffness 
will dramatically alter the low-temperature equilibrium 
structure of t-SHS polymers from that of the k^ = 
reference set, and that care should be taken in general- 
izing recent results obtained for flexible model polymers 
{e.g. Refs. [15-21]) to realistic polymers. In future stud- 
ies, we will enumerate and characterize the structural 
and mechanical properties of minimal energy semiflexi- 



Appendix A: Structure solver for packings with 
finite kc and kb 



We find the mechanically stable, minimal energy con- 
figuration {frnin} closcst to cach flexible packing {ffiex} 
by starting from configurations {fpert} representing small 
random perturbations [22] away from {rfiex}^ and mini- 
mizing t/chain- Energy minimization is carried out using 
the "dog- leg" method [23| , which is a combination of gra- 
dient descent and Newton's method [2^. Our numerical 
algorithm is schematically depicted in Figure [6l 

For finite the first spatial derivatives of t/chain 
evaluated at {r} ^ {rmin} are generally nonzero. We 
approach {r} = {fmin} using gradient descent and a 
second-order expansion of f{{r}) = t^chain({^})- This 
descent iteratively reduces / using 



fn= fn- (y^^fn, 



(Al) 



where fn is the value of / at the nth minimization step, 
/* is t/chain at the "knee" of the dog- leg [^] , the optimal 
step length is [26| 



a 



VfH{f){VfV' 



(A2) 



and H{fn) is the 3N x 3N Hessian matrix of fn- After the 
gradient descent step, we perform a least-squares Newton 
iteration on V/n using step- length control [24.] (Box 1): 

If the Frobenius norm ||V/n+i||2 of V/n+i decreases {i.e. 
||V/n+i||2 < ||V/*||2), we check whether /n+i < /*. If 
/n+i > /*, we set /n+i = /*; if /* > /n, the gra- 
dient descent has failed and we set /* = fn in Eq. 
[A3l We then check whether ||V/n+i||2 < ||V/*||2; if 
so. Equations lAltiASI represent one complete minimiza- 
tion step taking fn /n+i (Boxes 1-2). Otherwise, if 
||V/n+i||2 > ||V/*||2 the minimization algorithm has 
failed (Box 3); however, this occurs only for /c^ much 
larger than those considered here. 

The minimization procedure reduces the first spatial 
derivatives of t/chain to zero. When this condition is sat- 
isfied (V/ = 0), we check whether {r} represents a me- 
chanically stable solution, i.e. whether all 3N—6 nontriv- 
ial eigenvalues of H are positive (Box 4). Saddle-point 
solutions are avoided (when encountered) by perturbing 
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FIG. 6: Flowchart depicting the algorithm used to find the mechanically- stable, minimum-energy configuration {fmin} closest 
to each flexible packing configuration {ffiex}- 



coordinates along a negative eigenvector of the Hessian 
matrix (Box 5). 

Finally, should the solver reach the maximum num- 
ber of iterations for a particular {fpert} (MAXRUNS; 
Box 6) we attempt to minimize a different {fpert} ^ for 
up to MAXITS different {fpert} (Box 7). We found that 



MAXITS = 1000 different {fpert} with a maximum of 
50 iterations per {fpert} to be sufficient to minimize the 
energy or determine that no minimum that preserves the 
adjacency matrix exists for the N < 9 polymer packings 
considered here. 
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